Exactly solvable model of the 2D electrical double layer 
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We consider equilibrium statistical mechanics of a simplified model for the ideal conductor 
electrode in an interface contact with a classical semi-infinite electrolyte, modeled by the two- 
dimensional Coulomb gas of pointlike ± unit charges in the stability-against-coUapse regime of 
reduced inverse temperatures < [3 < 2. If there is a potential difference between the bulk interior 
of the electrolyte and the grounded interface, the electrolyte region close to the interface (known 
as the electrical double layer) carries some nonzero surface charge density. The model is mappable 
onto an integrable semi-infinite sine-Gordon theory with Dirichlet boundary conditions. The exact 
form-factor and boundary state information gained from the mapping provide asymptotic forms of 
the charge and number density profiles of electrolyte particles at large distances from the inter- 
face. The result for the asymptotic behavior of the induced electric potential, related to the charge 
density via the Poisson equation, confirms the validity of the concept of renormalized charge and 
the corresponding saturation hypothesis. It is documented on the non-perturbative result for the 
asymptotic density profile at a strictly nonzero /3 that the Debye-Hiickel /3 — > limit is a delicate 
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I. INTRODUCTION 

Asymmetric classical Coulomb mixtures, such as high- 
ly charged colloidal or polyelectrolyte suspensions, have 
attracted much attention in the last years due to the 
appearance of various anomalous phenomena; for a nice 
review, see Ref. m The contemporary theoretical treat- 
ment of equilibrium statistical mechanics of asymmetric 
Coulomb mixtures is based on the concept of renormal- 
ized charge. This concept has been introduced within the 
Wigner-Seitz cell models to describe an effective interac- 
tion between highly-charged "macro-ions" as a result of 
their strong positional correlations with the oppositely 
charged "micro-ions" HSk, 5]. The idea of renormal- 
ized charge is usually documented in the infinite dilu- 
tion limit on a simplified model of a unique charged col- 
loidal "guest" particle immersed in a symmetric weakly- 
coupled electrolyte 0, Q, 0, 0| • Plausible arguments were 
given to conjecture that the induced electric potential far 
from the guest colloid exhibits the form predicted by the 
Debye-Hiickel (DH) theory [sometimes called the linear 
Poisson-Boltzmann (PB) theory], with a renormalized- 
charge prefactor. Within the framework of the nonlinear 
PB theory, the renormalized charge saturates at some fi- 
nite value when the colloidal bare charge goes to infinity 

13 ■ A more general phenomenon of the saturation of 
the induced electric potential at each point of the elec- 
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trolyte region was studied in Ref. |Toi 

The idea of renormalized charge was developed within 
the linear and nonlinear versions of the PB theory, which 
are r igor ously valid in the limit of the infinite tempera- 
ture [ij. In order to treat correctly the Coulomb sys- 
tem at a finite temperature, one has to go beyond these 
mean-field approaches by incorporating electrostatic cor- 
relations among electrolyte particles. As concerns the 
large-distance asymptotic behavior of the electric poten- 
tial induced by a guest charge, this should affect both the 
renormalized-charge prefactor (provided that the concept 
of renormalized charge remains valid) as well as the corre- 
lation length of electrolyte particles which is expected to 
govern the exponential decay of the electric potential. In 
spite of the existence of many phcnomenological approx- 
imations based on heuristic extensions of the mean-field 
theories 0], there is no chance to solve exactly a three- 
dimensional (3D) Coulomb fluid at some finite tempera- 
ture. 

The situation is more optimistic in the case of 2D 
Coulomb fluids consisting of charged constituents with 
logarithmic pairwise interactions. These systems main- 
tain many generic properties, like screening and the re- 
lated sum rules [13, of "real" 3D Coulomb fluids. The 
2D Coulomb gas of symmetric ± unit point-like charges 
is stable against the collapse of positive-negative pairs of 
charges at high enough temperatures, namely for (3 < 2 
where (3 is the (dimensionless) inverse temperature or 
the coupling constant. The collapse starts to occur at 
j3 = 2: interestingly, although the free energy and each 
of the species densities diverge, the truncated Ursell cor- 
relation functions are finite at this inverse temperature. 
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The collapse /? = 2 point is exactly solvable due to its 
equivalence with the free-fermion point of the Thirring 
fermionic representation of the 2D Coulomb gas [T^ IT3 | . 
The exact solution involves the bulk thermodynamics for 
an infinite system and special cases of the surface ther- 
modynamics for a semi-infinite Coulomb gas in contact 
with an impermeable dielectric wall; for an exhaustive 
review, see Ref. [H 

Besides the Thirring fermion representation, the 2D 
Coulomb gas is equivalent to the 2D Euclidean, or (1+1)- 
dimensional quantum, sine-Gordon field theory with a 
conformal normalization of the cos- field [3. Although 
the bulk 2D sine-Gordon model has been known to pos- 
sess the integrability property from the seventies |l7j | , the 
explicit exact solution of its ground-state characteristics 
was derived only quite recently due to a progress in the 
method of Thermodynamic Bethe Ansatz (TBA) |Tl[l9l|. 
Based on the equivalence with the sine-Gordon theory, 
the bulk thermodynamic properties (free energy, internal 
energy, specific heat, etc.) of the 2D Coulomb gas have 
been obtained exactly in the whole stability region of 
point-like charges (3 <2 [23| ■ From a gnoseological point 
of view, this is the only exactly solvable case of a contin- 
uous (i.e. not on a lattice) fluid in more than one dimen- 
sion. Later on, the form- factor approach ^21.. .22.] was ap- 
plied to calculate the large-distance asymptotic behavior 
of the charge and number density |3| pa-ir correla- 
tion functions in the bulk of the infinite 2D plasma. The 
integrability of the 2D sine-Gordon theory was shown to 
be preserved also in the half-space geometry with specific 
types of boundary conditions jl^. The surface thermo- 
dynamic properties (surface tension) of the 2D Coulomb 
gas in contact with an ideal conductor wall 26] and an 
ideal dielectric wall [23| were obtained in the whole sta- 
bility range of /3 < 2 through the mappings onto the 
boundary sine-Gordon model with Dirichlet and Neu- 
mann boundary conditions, respectively, with the aid 
of the known TBA results |2^ for these integrable 
boundary field theories. 

The exact treatment of the 2D Coulomb gas was re- 
cently extended also to inhomogeneous situations of the 
present interest when one guest arbitrarily-charged par- 
ticle is immersed in the bulk of an electrolyte modeled by 
the 2D Coulomb gas. At the free-fermion (collapse) point 
(3 — 2, the considered problem is solvable in the Thirring 
format ^3 even for the guest Q-charged particle being 
of colloidal type, i.e. possessing a hard core of radius a 
which is impenetrable to the electrolyte ± unit charges. 
Based on an explicit formula for the electric potential 
induced by the charged colloid in the electrolyte region, 
the concept of the renormalized charge 0, S S was 
shown to fail in this strong-coupling regime. On the other 
hand, in the limit Q ^ 00, the anticipated phenomenon 
of the electric potential saturation [l^ was confirmed at 
this free-fermion point. The special case of the point-like 
(a = 0) guest Q-charge is solvable, within the framework 
of the sine-Gordon format, inside the whole stability in- 
terval of the electrolyte (3 < 2 Isil]. The explicit results 



for the asymptotic behavior of the induced electric poten- 
tial confirm the adequacy of the concept of renormalized 
charge in this weak-coupling regime. The exact results 
are rigorously valid provided that (3\Q\ < 2, i.e. when the 
guest Q-charge docs not collapse with an opposite unit 
charge (counterion) from the electrolyte. The possibil- 
ity of an analytic continuation of the results beyond the 
stability border f3\Q\ = 2 was conjectured 31], however, 
the validity of this "regularization hypothesis" was not 
rigorously proved. The restricted rigorous validity of the 
exact results to the region f3\Q\ < 2 prevents one from 
studying the saturation phenomena in the limit Q — > 00. 

In order to avoid the collapse of point-like guest charges 
with electrolyte counterions, one has to search for an- 
other integrable 2D model with the guest charges uni- 
formly smeared over a line manifold. The simplest system 
of this kind is the 2D half-space Coulomb gas in contact 
with a plain hard wall carrying a uniform "line" charge 
density. Although this model is exactly solvable at the 
free-fermion (3 — 2 point ^3] j it can be easily shown that 
its sine-Gordon formulation does not belong to the fam- 
ily of the boundary sine-Gordon theories integrable at 
arbitrary (3 < 2 |25j. Another way how to introduce an 
interface charge density is to consider a simplified model 
for an electrode in contact with a classical electrolyte: 
the half-space Coulomb gas bounded by an ideal con- 
ductor wall, with a potential difference if between the 
bulk interior of the electrolyte and the grounded interface 
[3^ Is^ . As soon as (/3 ^ 0, the region of the Coulomb 
fiuid close to the interface (known as the electrical double 
layer) carries some nonzero surface charge density. The 
2D version of the model was mapped in Rcf. "2? onto the 
integrable half-space sine-Gordon model with a specific 
(/^-dependent Dirichlet boundary condition. 

The present paper concentrates just on this integrable 
model of the 2D electrical double layer and has two main 
aims. The first aim is rather technical: we present a 
method for obtaining the charge and number density pro- 
files of electrolyte particles at asymptotically large dis- 
tances from model's interface. This task is equivalent 
to the calculation of one-point functions of bulk fields in 
the boundary sine-Gordon model with Dirichlet bound- 
ary conditions. Although the one-point functions of cer- 
tain specific boundary field theories have already been 
analyzed using the truncated conformal space approach 
and the form- factor expansion [ssl l36| , some additional 
generalizations and technicalities have to be developed 
for the present model. Among others we document on 
the exact non-perturbative asymptotic form of the num- 
ber density profile at large distances from the interface 
that the DH /3 — > limit is a delicate point which has to 
be taken with cautiousness. The second aim consists in 
pointing out physical consequences of the obtained exact 
results. The knowledge of the asymptotic behavior of the 
charge density at large distances from the interface en- 
ables us to derive the asymptotic large-distance tendency 
of the induced electric potential toward its bulk value Lp. 
In the whole stability range of the electrolyte coupling 
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P < 2, the asymptotic form of the electric potential co- 
incides, up to a renormalized-charge prefactor and the 
plasma correlation length, with the one obtained in the 
DH limit. This result supports the general validity of 
the concept of renormalized charge in the weak-coupling 
regime of the electrolyte. The saturation hypothesis of 
the induced electric potential is also confirmed. 

The article is organized as follows. In Sec. II, we intro- 
duce the notation and briefly summarize important as- 
pects of the mapping of the infinite 2D Coulomb gas onto 
the bulk (1+1 )-dimensional sine-Gordon theory. Sec. Ill 
deals with the 2D electrical double layer of interest and its 
mapping onto the semi-infinite sine-Gordon theory with 
Dirichlet boundary conditions. We would like to empha- 
size that the presentation in Sees. II and III is sketchy; 
for a detailed study of some specific points, we quote rele- 
vant references. Sec. IV presents the mean-field theories 
for the 2D electrical double layer: the DH /3 ^ limit 
with the leading /3-correction, and the nonlinear PB the- 
ory. The crucial Sec. V is devoted to the derivation of 
the asymptotic charge and number density profiles for the 
2D electrical double layer, in the whole stability region of 
inverse temperatures /? < 2. The Thirring free-fermion 
point /3 = 2 is discussed in Sec. VI. A brief recapitulation 
and some concluding remarks are given in Sec. VII. 



II. THE 2D COULOMB GAS IN THE BULK 

A. Basic definitions 

We start with a brief description of the classical 
Coulomb gas formulated in the infinite 2D space of points 
r G M^. It is realized as the limit of a finite system with 
periodic boundary conditions. The system consists of 
point-like particles {i} of charge {qi = ±1} (the elemen- 
tary charge e is set for simplicity to unity) immersed in a 
homogeneous medium of dielectric constant e = 1. The 
interaction energy E of a, set of particles {qi, r^} is given 
by 



E{{q^, r J) = ^ qrqjv{\r, ~ r 



(2.1) 



where the electrostatic potential v is the solution of the 
2D Poisson equation 



Av(r) = -27r(5(r) 

subject to the boundary condition Vw(r) 
oo. Explicitly, one has 



r e 



as 



v{r) = - In ( — - ) , 



(2.2) 



(2.3) 



The free length constant tq, which fixes the zero point 
of the Coulomb potential, will be set for simplicity to 
unity. The Fourier transform of the 2D Coulomb po- 
tential (|2.3|l exhibits the form l/|kp with the charac- 
teristic singularity at k — > 0. This maintains many 



generic properties (like screening and the related sum 
rules [IJ) of 3D Coulomb fluids with the interaction po- 
tential ?;(r) = l/|r|, r e M^. 

The system is treated in thermodynamic equilibrium, 
via the grand canonical ensemble characterized by the 
(dimensionless) inverse temperature /3 and the couple of 
particle fugacities z+ = z_ = z. The grand partition 
function is defined by 



N+=0 N-=0 



N+l iV_ 



(2.4a) 



where 



N 



Q{N+,N_)= J] d^r.exp [-/?£({(?„ r J)] (2.4b) 



is the configuration integral of iV_|_ positive and nega- 
tive charges, and N = N+ + For the considered case 
of point-like particles, the singularity of the Coulomb po- 
tential (|2.3|l at the origin r = can cause the thermody- 
namic collapse of positive-negative pairs of charges. The 
stability regime against this collapse is associated with 
the 2D spatial integrability of the corresponding Boltz- 
mann factor exp[/3ti(r)] = |r|~'^ at short distances, and 
therefore corresponds to small enough inverse tempera- 
tures (i < 2. At large distances, the Coulomb interac- 
tion is screened to a short-distance effective interaction 
of Yukawa type. The infinite system is homogeneous and 
translationally invariant. Denoting by (• • • )^ the thermal 
average, the number density of particles of one charge 
sign q (= ±1) is defined by Uq = " i"i))/3- 



Due to the charge symmetry, 



n/2 where 



n is the total number density of particles. At the 
two-particle level, one introduces the two-body densities 
"?g'(|r-i"'|) = (Ei5:^j^M,^(i"-r»)V,9,^(r'-i-j))/3, etc. 



B. The sine-Gordon representation 

The infinite 2D Coulomb gas is mappable onto the bulk 
sine-Gordon theory Using the fact that, according 
to Eq. (|2.2|l . — A/(27r) is the inverse operator of the 
Coulomb potential v and renormalizing the particle fu- 
gacity z by the (divergent) self-energy term exp[/3'i;(0)/2], 
the grand partition function (|2.4() can be turned via the 
Hubbard-Stratonovich transformation into 



S(z) 



/ X'0exp[-S'(z)] 
/X>0exp[-S'(O)]' 



(2.5a) 



where 



S{z) 



dV 



-^(V0)'-2zcos(60)l , 6^ = ^ 
Idtt J 4 

(2.5b) 

is the Euclidean action of the 2D sine-Gordon model. 
Here, 0(r) is a real scalar field and / 1?0 denotes the 
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functional integration over this field. The one- and two- 
particle densities are expressible as averages over the sine- 
Gordon action (|2.5b|l in the following way 



Defining the "topological charge" q as 



'(|r-r'|) 



z,(e'«''*), 



(2.6a) 
(2.6b) 



The renormalized fugacity parameter z gets a precise 
meaning when one fixes the normalization of the cou- 
pled cos-field. In the Coulomb-gas format, the short- 
distance behavior of the two-body density for the op- 
positely charged particles is dominated by the corre- 
sponding Boltzmann factor of the Coulomb potential: 

(r, r') ~ z+Z-\y — y'\~^ as |r — r'| 0. With respect 

to the representation l|2.6b|l . in the sine-Gordon picture 
this corresponds to 



r — r' 



as r — r 



0. (2.7) 



Under this short-distance (conformal) normalization, the 
divergent self-energy factor disappears from statistical re- 
lations calculated in the sine-Gordon format. 

The 2D Euclidean sine-Gordon action (|2.5b|l takes its 
minimum at a (/)(r) constant in space. Due to a dis- 
crete symmetry — > (/> -|- 27rri,/6 (n being an integer), 
the action has infinitely many ground states |0„) char- 
acterized by the associate expectation values of the field 
((/))„ = 27rn/6. In the considered infinite-volume limit 
and for &^ < 1, these ground states become all degener- 
ate 0. Equivalently, the discrete (/)-symmetry is spon- 
taneously broken an so it is sufficient to develop the sine- 
Gordon action Ij2.5b|) around any one of its ground states, 
say the one |0o) with ((/))o = 0. 

In order to pass from the present Lagrangian formu- 
lation (|2.5|) to the Hamiltonian one, one chooses say the 
x-direction to be the "Euclidean time" and associates a 
Hilbert space TL to any equal-time section {x = const, y e 
(— cx),cx))}. (Choosing the y-direction to be the "Eu- 
clidean time" gives equivalent quantization). General 
states are vectors in TL whose evolution is governed by 
the Hamiltonian operator 



H = 



dy 



47rn2 + —(dycpf + 2zcos(b(f>) 
Ion 



(2.8) 



In the considered region < 1 with the spontaneously 
broken discrete ^-symmetry, the sine-Gordon field the- 
ory is massive in the sense that Ti, is the Fock space of 
multiparticle states. After rotation x — it to the (1-1-1) 
Minkowski time-space {t,y), these multiparticle states 
are interpreted as the asymptotic "in-" and "out-" scat- 
tering states (see below). 

The particle spectrum of the ( 1 + 1 )-dimensional sine- 
Gordon field theory is the following [13. The basic par- 
ticles are the soliton S and the antisoliton S which form 
a particle-antiparticle pair of equal masses M. They cor- 
respond to specific (^-configurations that interpolate be- 
tween two neighbouring ground states, say |0o) and |0i). 



b 

2^ 



<^y-T-(t>{x,y) 

, oy 



b 

2^ 



{x, co) — 4>{x, —00)] 



(2.9) 

g = +1 (—1) for the soliton (antisoliton). Due to topo- 
logical reasons, the soliton and the antisoliton can coexist 
in the particle spectrum only in neutral pairs. The S — S 
pair can create neutral {q = 0) bound states {Bi]i = 
1,2,... < p^^}; these particles are called "breathers". 
Their number depends on the inverse of the parameter 



1-52 



4-/3 



The mass of the Bi-breather is given by 



nil 



2M sin 



(2.10) 



(2.11) 



and this breather disappears from the particle spectrum 
just when = 2M (i.e. p — 1/?). Note that the 
breathers exist only in the stability region of the point- 
like Coulomb gas < /3 < 2 (0 < p < 1); the fightest Bi- 
breather disappears just at the border 6^ = 1/2 (/3 = 2), 
which is the field-theoretical manifestation of the collapse 
phenomenon. The S — S pair remains in the spectrum up 
to the Kosterlitz-Thouless transition point b^ = 1 {P = A) 
beyond which the sine-Gordon model ceases to be mas- 
sive. In what follows, we shall restrict ourselves to the 
stability region of point- like charges (3 < 2. 

Let a e {5, S, Bi{i = 1, 2, . . . < p~^)} denote the type 
of the given particle and the corresponding particle 
mass. Since the sine-Gordon model is a relativistic field 
theory, the energy E and the momentum p of the particle 
can be parameterized as follows 



Ea = rua cosh0. Pa — TUa sinh0. 



(2.12) 



where G (—00,00) is the particle rapidity. The asymp- 
totic n-particle states { |n) } are generated by the "particle 
creation operators" A'^{9), 

<(^?iX(^2)---A+ ((?„)|0), (2.13) 

where |0) G is the ground state of the Hamiltonian 
H given by H2.8|l . The state H2.13|l is interpreted as an 
in-state if the rapidities are ordered as 6*1 > 6*2 > • • • > 
On and as an out-state in the case of the reverse order 
9i < 62 < ■ ■ ■ < On- The in- state basis and the out-state 
basis are related via the scattering n ^ n S'-matrix. The 
2 — !■ 2 process is described simply by 

At^ieMtM) - E s'atlio^ - ^2)A+(^^2)A+ (0i). 

(2.14) 

Here, the momentum conservation demands rriai = "^61 



and nia 



, so that the inequalities ai 7^ 61 or 



0-2 7^ ^2 are allowed only for the soliton-antisoliton pair 
with the degenerate masses equal to M. Like for any 
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integrable field theory, the n ^ n S'-matrix of the sine- 
Gordon model factorizes into a product of n{n — l)/2 
two-particle 5'-matrices. The two-particle S'-matrix pos- 
sesses many symmetry constraints and its explicit form 
was obtained by exploring four general axioms: the Yang- 
Baxter equation; a unitarity condition; analyticity and 
crossing symmetry; the bootstrap principle MT|. 

The (dimensionless) specific grand potential uj of the 
sine-Gordon model H2.5() . defined in the infinite- volume 
limit as 



1 



InS, 



was found by using the TBA in Ref. |T| 



8 sin(7rp) 



(2.15) 



(2.16) 



Here, mi is the mass of the lightest J5i-breather [see for- 
mula 12.11|l taken with i = 1] and the parameter p is 
defined by (|2.in(l . Under the conformal normalization 
H2.7|l . the relationship between the fu gaci ty z and the 
soliton mass M was established in Ref. lift 



r(62 



7rr(l - fe2) 



M 



V^r((i + p)/2) 
2r(p/2) 



2-26^ 



(2.17) 



where F stands for the Gamma function. Since the total 
particle number density n of the 2D Coulomb gas is given 
by the obvious equality 



dz ' 



(2.18) 



Eqs. H2.16|l and H2.17|l imply the explicit density- fugacity 
relationship, and consequently the complete bulk thermo- 
dynamics, of the 2D Coulomb gas in the whole stability 
region (3 < 2 The mass mi plays the role of the in- 
verse correlation length of the plasma particles [2^ |2^ . 
Using Eqs. (|2.16|l - (|2.18|l . it is expressible as 



mi 



sin(7rp) 



-,1/2 



192 



where 



K = y/2Tr(3r 



is the inverse Debye length. In the DH limit /? 
reduces to k as it should be. 



(2.19) 

(2.20) 

0, mi 



r € defined by Cartesian coordinates (x,y). The 
electrode-electrolyte interface is localized at a; = 0, along 
the y axis. The half-space A — {{x,y);x < 0} is oc- 
cupied by an ideal-conductor wall of dielectric constant 
ew — ^ oo, impenetrable to electrolyte particles. The 
electrolyte is localized in the complementary half-space 
A = {{x,y);x > 0}. It is modeled by the classical 2D 
Coulomb gas of point-like unit ± charges. The interface 
{x — 0) is kept at zero potential while the bulk interior of 
the electrolyte (x ^ oo) is assumed to be at some elec- 
trostatic potential (p. The non-zero potential f causes a 
splitting of the charge fugacities: 



z± — z exp (±/3(p) . 



(3.1) 



Alternatively, chemical potentials /i+ and /i_ can be de- 
fined by z± = exp{(3^±)/X'^ where A is the de Broglie 
thermal wavelength. The bulk Coulomb gas is neutral 
|37l |. and therefore its bulk properties depend only on 
the chemical potential combination fi = (/i_|_ -I- /Lt_)/2, 
i.e. on z. The difference /i+ — /i_ (or (p) is relevant only 
for the surface properties of the electrolyte region close 
to the interface (the electrical double layer) [s^. l3^ ; if 
tp ^ the electrical double layer carries some surface 
charge density. 

The presence of the ideal-conductor wall is described 
mathematically by charge images: the particle with 
charge q localized in the electrolyte region at the point 
r = (x > 0,y) induces the image with the opposite 
charge q* = —q localized in the wall region at the point 
r* = {—X, y). The interaction energy E oi a. set of parti- 
cles {qi,ri — {xi > 0, yi)} then consists of two parts (see, 
for example, Ref. ISa) : 

E{{q^,Yi}) ^^q,qjv{\vi - Yj\) + ]^'^q,q*v{\vi ~ r*|); 

the first term corresponds to direct particle-particle in- 
teractions, while the second term describes interactions 
of particles with the images of other particles and with 
their self-images. The grand partition function Sbry is 
again given by 



oo oo Af+ N- 



Ar^=0 Af_=0 



N+\ NJ. 



(3.3a) 



where the configuration integral 



TV 



Q(iV+,iV_)= / []d2r,exph/3S({<7„rJ)] (3.3b) 



III. THE 2D ELECTRICAL DOUBLE LAYER 

A. The definition 

We first introduce in detail the model of interest which 
describes an electrode in contact with a classical elec- 
trolyte. Let us consider an infinite 2D space of points 



is now restricted to the half-space x > Q. The stability 
range of inverse temperatures for the surface thermody- 
namics is determined by the Coulomb interaction of the 
charged particle with its image |0|: the Boltzmann fac- 
tor of a particle at a distance x from the wall with its own 
image, proportional to {2x)~^/'^^ is integrable at small ID 
distances x provided that (3 <2. Note that the bulk and 
surface thermodynamic stability intervals of {3 coincide. 
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Due to the translational invariance of the system along 
the y axis, the species densities depend only on the a;- 
coordinate: n±(r) = n±{x), x > 0. The electrolyte 
is neutral in the bulk interior of the electrolyte, i.e. it 
holds \vaix~too n±{x) — n± = n/2. The species densities 
determine the particle number density 

n{x) = n,q{x) — n+{x) +n-{x) (3-4) 

q=±l 

and the charge density 

p{x) — qnq{x) = n+{x) ~ n-{x). (3-5) 

q=±l 

The induced (averaged) electrostatic potential (p{x) in 
the electrolyte region is related to the charge density via 
the Poisson equation 



(p"{x) = -27rp(x). 



(3.6) 



The potential satisfies the obvious boundary conditions 
(p{0) = and (p{oo) = 95, and the regularity requirement 
(all derivatives if' , ip" , etc. vanish) at a; — > 00. Since we 
will be interested in the asymptotic approach of (p{x) to 
its bulk value ip at large distance x from the electrode 
surface, it is natural to introduce the quantity 



S(p{x) — (p{x) — ip 



(3.7) 



which vanishes as x ^ 00. 

There are two sum rules which can be derived without 
solving explicitly the boundary problem. First, the con- 
sideration of Eq. (|3.t)|) in the integral dx p{x) implies 



dxp(x) = — vp'(O). (3.. 



We note that a nonzero surface (more precisely "line") 
charge in the electrolyte /p dxp(x) can appear only in 
the special case of an ideal-conductor wall where this 
charge is exactly compensated by the opposite surface 
charge a of particle images, and the system as a whole is 
neutral. Since an isolated charged line at x = 0, carrying 
a uniform charge a per unit length, induces the electric 
potential such that 



^'(0) = -27rcr, 



(3.9) 



the sum rule (|3.8|) can be understood as the neutrality- 
type condition 



/•OC 

/ dx p{x) + o- = 0. 
Jo 



(3.10) 



Second, the consideration of Eq. H3.6|l in the inte- 
gral /p Ax xp{x) and the subsequent two integrations by 
parts lead to 



j dxxp{x)^—, 



(3.11) 



i.e. the dipole moment of the charge density is related 
to the potential difference across the electrical double 
layer. This relation follows from elementary macroscopic 
electrostatics 38] . 

At small distance from the wall a; 0, the species 
densities are determined by the Boltzmann factors of the 
corresponding particle with its image (2^: 



n±{x) 



z± 



{2xY/^ ' 



0. 



Consequently, 



cosh(/3(/?) 
n[x) ^ 2z— — ttttt:-. 



(3.12) 

(3.13a) 
(3.13b) 



B. The boundary sine-Gordon representation 

The considered particle- image system is mappable onto 
a boundary sine-Gordon theory |26| |. In particular, the 
grand partition function H3.3|l can be written as 



J Vcj) exp (-5'bry(^:)) 
/ I?0exp(-5bry(O))^ 



(3.14a) 



where 



d^r 



a:>0 



-^(V(/.)'-2zcos(6</)) 



(3.14b) 

is the 2D Euclidean action of the boundary sine-Gordon 
model defined in the half-space x > and the real scalar 
field (/)(r) fulfills the following Dirichlet conditions at the 
X — Q boundary: 



(j){x = 0,y) = (/lo 



-4i6i^. 



(3.15) 



The fact that the boundary value of the field, (/>o , is a pure 
imaginary number makes no problem: as is usual in the 
field theory, first one expresses the quantities of interest 
as functions of real 0o a-nd then analytically continues the 
obtained results to complex values of (pQ. This procedure 
was successfully applied to the calculation of the surface 
tension (i.e. the surface part of the grand potential) for 
the present model j2^]. 4>q will usually appear in the 
combination 



= 2i^. 
2b ^ 



(3.16) 



Within the formalism developed in Ref. l26l the one- 
particle densities n± (x) in the electrolyte region x > are 
expressible as averages over the boundary sine-Gordon 
action (|3.14b|) with Dirichlet boundary conditions (|3.15|) 
as follows 



n±(x) 



z(e±'''^(-.«))bry. 



(3.17) 
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Here, regarding the y-invariance of the boundary mean 
values (e''='''"*^''^)bry, we have set y = for simpUcity. 
Thus, the particle number density H3.4|l and the charge 
density (|3.5|l read 



)bry (*5 



-if)0(2:,O)\ 



bry 



, (3.18) 
, (3.19) 



respectively. 

The specific case of Dirichlet boundary conditions in 
the semi-infinite sine-Gordon model does not spoil the 
integrability property of the bulk theory 25] . In passing 
from the Lagrangian formulation H3.14|l to a Hamilto- 
nian one, contrary to the bulk case, we have two different 
choices considering either a; or y as the Euclidian time. 

If y is taken to be the Euclidean time, then the bound- 
ary is in space, and a boundary Hilbert space Ti^ is as- 
sociated to any time slices {x G {Q,oo),y ~ const}. The 
time evolution is governed by the Hamiltonian 



Hb = 



dy 



2zcos(6(/)) 



(3.20) 

where the boundary condition, (|3.15|l . is satisfied. The 
boundary Hilbert space consists of boundary bound 
states, |34| . and bulk multiparticle states: 

<(0i)<(e2)---< (ej|0)B, (3.21) 

where |0)b eHb is the ground state of the Hamiltonian 
Hb given by H3.20(l . The state H3.21|l is interpreted as an 
in-state if the rapidities are ordered as > > 6*2 > 
• • • > 9n and as an out-state in the case of the reverse 
order < 6i < 62 < • ■ ■ < On- The in-state basis and 
the out-state basis are related via the reflection n — > 
n i?-matrix. This matrix factorizes into the product of 
pairwise scatterings H2.14II and of, R'^{9)-s, the individual 
one-particle reflections A'^{—9) — > ^h^(^) {9 positive) off 
the boundary. Owing to the energy conservation, Ra{9) 
vanishes if ^ mi,. The i?-amplitudes were obtained 
explicitly for the soliton/antisoliton pair in Ref. 12^ and 
for the breathers in Ref. We shall need the Dirichlet 
i?-amplitudes for the lowest Bi- and B2-breathers; with 
the notation Rg^_{9) = Rg\9), one has explicitly 

l)(f + l)(f-^) 



and 



R^b\0) 



(3.22a) 



;i)(i-fl)(p+l) 



+ 



where we used the symbol 

sinh(f-K^) 



ix) 



sinh ( 



U_ 17TX ) 

2 2 J 



(3.22b) 



(3.23) 



The reflection amplitudes of breathers have simple poles 
at the imaginary rapidity 9 = m/2. In the particular 
case of the first two breathers H3.22|l . one finds that 



R^\9) 



29- 



(3.24) 



where the "boundary couplings" gi and 52 are extracted 
in the form 



gi — 2 tan 



1 



l_cos(If)+sin(2f) 



1/2 



92 



2tan(2f 



HI) tan (if 



£21 1 

2 J 



tan(^) [tan(if)tan(f + 2f)] 



1/2- 



, (3.25a) 
(3.25b) 



In the alternative Hamiltonian description one can take 
X to be the Euclidean time and associate with any equal- 
time section {x = const, y € (— oo,cxd)} the same Hilbert 
space Ti. as in the bulk theory. The Hamiltonian operator 
is now given by Eq. H2.8|l . The boundary at x = ap- 
pears as the initial condition described by the boundary 
state \B) G Ti.. In the 1-1-1 Minkowski space-time, \B) 
can be written as a superposition of the bulk asymptotic 
states (gusli: 



\B) = expJ ^5,AJ(0) 

I a 



d9_ 
2tt 



Y,K^\0)A+i-9)Ati9)\\0). (3.26) 



(a.b) 



The amplitude K°-^{9) is related to the reflection matrix 



K''\9) = Rl 



(3.27) 



where a denotes the antiparticle of the particle a {Bj — 
Bj). Ghoshal and Zamolodchikov [2^ identified cja with 
the boundary coupling ga, but Dorey et al 35] found in 
the case of the Lee- Yang model the relation 



Qa = gafi- 



(3.28) 



The strong evidence that this formula extends to any 
boundary (l-l-l)-dimensional quantum field theory, and 
in particular to the sine-Gordon one, was presented later 
in Ref. ,3^ 



IV. MEAN-FIELD THEORIES 
A. The Debye-Hiickel limit 

The position-dependent species densities in the elec- 
trical double layer can be evaluated systematically via 
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an expansion in powers of /3 around the DH high- 
temperature hmit /3 — *■ 0. The technique is based 
on a Mayer expansion of the free energy with series- 
renormahzed bonds between each couple of field circles; 
for a detailed description of the method see Ref. 

In the lowest expansion order, the species are consid- 
ered to be constant in the electrolyte region, n±(x) — n/2 
{x > 0). Under this assumption, the renormalized bond 
IS given by if = X(o) with 

K^''Hrur2) = -fiKoi^in^) + PKoinrl^). (4.1) 

Here, k is the inverse Debye length (|2.20() . Kq is the 
modified Bessel function of second kind, ri2 = |ri — 
and = ~ — ~ ^2\- The first correction to 
the constant species densities can be obtained iteratively 
by inserting the lowest-order if^"^ into a basic generating 
formula for the density profiles, with the result 

n±{x) = - exp <^ ±l3[ip - ip{x)] -f -Ko{2h 



Ikx) \ . (4.2) 



Here, f{x) is the electric potential induced by the charge 
density p{x) via the Poisson equation H3.6|l . 

Expanding the exponential in H4.2|) to order (3 gives for 
the charge density 



B. The leading high-temperature correction 

The next iteration for n±{x) is obtained by using the 
DH density l|4.7(l in the basic equation for the renormal- 
ized bond K and treating 6n{x) as a perturbation. Now 
K = iff") + if(i), where iff"' cx /3 is defined by Eq. 
(14. 1|) and K^^^ cx is given, to first order in the density 
perturbation, by the integral equation 

if(i)(ri,r2)= / d\3K^°Hrur3)6n{xs)K'^"\r,,T2). 

Jx3>0 

(4.9) 

Taking also K'-^^ into account, instead of 14.2|l we have 



TP5ip{x) + ^Ko{2Kx) + ^K('\x)^ 



p{x) ^ (3n[ip- ip{x)] 



(4.3) 



(4.10) 

where i^T*^^' (r, r), with r having the coordinate x, is re- 
named K^^^{x). Expanding the exponential in (|4.1U|) up 
to order /3^, one obtains 

p{x) = -(3n6(p{x) - -^Koi2Kx)Sip{x), (4.11) 
6n{x) = I^Ko{2kx)+'^K^'Hx) 

+ ^[6^ix)f + ^[Koi2Kx)r. {4.12) 



Inserting this into Eq. H3.6I) and considering the bound- 
ary conditions for the electric potential, one gets 

dipDii{x) = -ipexp{-Kx), (4.4a) 
Punix) = f3nipe:icp{—Kx), (4.4b) 

where the deviation of the electric potential from its bulk 
value is defined by Eq. (|3.7fl . The surface DH "image- 
charge" (Tdh, defined by the couple of equivalent Eqs. 
(|3.9|) and (|3.1U|) . is obtained in the form 



fDH 



Kip 

2^' 



(4.5) 



In terms of ctdh , the relation 14.4al) is written as 



Sfuaix) = ^exp(— Kx). (4-6) 



Inserting the charge density H4.11|l into the Poisson 
equation (|3.6(l and considering the /3-cxpansion in the 
form 

S^{x) = -^e-«- + ^f(Kx), (4.13) 

the unknown /-function is determined by the ordinary 
differential equation 



fix) - fix) = -e--Ko{2x) 



(4.14) 



with the zero boundary conditions /(O) = /(oo) = 0. In 
terms of the ID Green function 



G{x,x') = -e ^>sinh(a;<) 



(4.15) 



with x< = min{a;, x'} and x> = max{a;, x'}, the solution 
of Eq. reads 



Expanding the exponential in (|4.2|l to order /?, the par- 
ticle number density reads 

(3n 

n-Dnix) = n+Snouix), Snuaix) ^ —Ko{2kx). (4.7) 
At large distances x from the wall, 

(5nDH(a;) ~ (^^^j exp(-2Ka;), x oo. (4.8) 



/(x)= / dx'G{x,x') -e-""' Kq{2x') 



(4.16) 



At large x, f{x) behaves like e ^[{tt/2) — l]/4. Conse- 
quently, 



5ip(x) 
p{x) 



-if 



(3nLp 



1-^f^-l 
8 V2 

1-^f^-l 
8 V2 



(4.17a) 
(4.17b) 
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These formulae provide the pure exponential asymptotic 
decay of the DH results H4.4I I. in agreement with the con- 
cept of renormalized charge 0, Q, H 0| ■ The renormalized 
image charge tTrcn, defined in analogy with Eq. (|4.6II as 



2-Kar, 



■ exp(— ra), 



has the following weak-coupling expansion 

11 +0(/3') 



'2^ 



8 V2 



(4.18) 



(4.19) 



From Eq. H4.16|l one gets that /'(O) = 1/2. Combining 
then Eqs. H3.9|l and (|4.13|) . the leading /3-correction to 
the surface "image-charge" ()4.5|l reads 



"2^ 



l + ^ + OiP') 



(4.20) 



The problem of the number density deviation from its 
bulk value, Sn(x) given by Eq. (|4.12l) . is more compli- 
cated because the function K'^^\x) is only defined im- 
plicitly as the solution of the integral equation H4.9|l . It 
can be shown after lengthy algebra that 



16 



exp(— 2Ka;) as x — > oo. 



Thus, at large distances from the wall, 

bn[x) ~ — - + T77 exp(-2Kx). 



(4.21) 



(4.22) 



Note that this asymptotic behavior differs fundamentally 
from, and is superior to, that obtained in the DH limit 
(I4.8|l . We conclude that, in contrast to the charge density, 
the DH theory does not provide an adequate description 
of the large-distance decay of the particle number density 
to its bulk value. 



C. The nonlinear Poisson-Boltzmann theory 

In the nonlinear PB theory, one keeps the deviation of 
the electrostatic potential from its bulk value, 1^9(2;) — lys, 
in the exponential form (|4.2() . The DH expression for the 
charge density 1)4. 3|) then takes the nonlinear form 

p(x) = nsinh{/?[(p - (^(x)]}. (4.23) 

The corresponding Poisson equation 

^J'{x) ^ -27rnsinh{/3 - Lp{x)\} (4.24) 

is subject to the obvious boundary conditions <y5(0) = 
and <^(oo) = kp. We recall the well-known fact that, with 
the identifications <^stix") = — 4i6(y5(a;) and z — the 
nonlinear PB equation 14.24|l corresponds to the static 
equation (0 does not depend on "time" y) of the "classi- 
cal" variational treatment of the boundary sine-Gordon 
action Ij3.14b|l . 5S\:,'cy[4'st{x)] = 0. 



For the present semi-infinite geometry, the solution of 
the nonlinear PB equation 14.24|l can be derived explic- 
itly: 



5ip{x) 



f3' 



■In 



l + c-'"^tanh(/3(/3/4) 



1 - e-'*^tanh(^(p/4) 
At asymptotically large distance x from the wall, 

S^{x) - -Itanhf 



(4.25) 



(4.26) 



This behavior is again in full agreement with the idea 
of renormalized charge. The renormalized image charge 
(Tien, defined by Eq. 1)4.18(1 . reads 



2k [pip 

Crcn = ^ tanh I — 

np 



(4.27) 



In the (3^0 limit and for a finite (p, the DH result H4.5() 
is reproduced as it should be. On the other hand, the 
leading /3-correction to the DH result in Eq. H4.19|l is not 
reproduced. We shall show in the next section that the 
nonlinear formula H4.27|l describes correctly the scaling 
regime of limits /3 ^ and ip ^ 00 with the product P(p 
being finite. 

The surface "image-charge" a, defined by either H3.9fl 
or H3.10|l . is obtained as follows 



7r/3 \ 2 



(4.28) 



For a nonzero /3 and in the limit of the infinite cr-charge, 
which is equivalent in view of Eq. 14.28|1 to the limit 
ip ~* 00, the renormalized a^cn (I4.27|) saturates at the 
finite value 



2k 



(4.29) 



More generally, for a nonzero (3 and in the limit cp 
Eq. (|4.25(l reduces to 



00, 



Sip*ix) 



■In 



1 + CXp( — KX) 

1 — exp(— Kx) 



(4.30) 



i.e. the electric potential deviation from its bulk value 
saturates at a finite value in every point of the electrolyte 
region a; > 0, in accordance with the hypothesis of the 
potential saturation 1^ . Note that the potential satura- 
tion is a pure nonlinear effect: there is no saturation in 
the DH relation H4.4a|l . 



V. ASYMPTOTIC CHARGE AND NUMBER 
DENSITY PROFILES FOR /3 < 2 

A. Boundary one-point functions: general 
formalism 

Let us consider a general integrable 2D boundary field 
theory, defined in the half-space a; > and possessing an 
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integrable boundary condition at a; = 0. For the time 
being, the spectrum of the corresponding bulk theory is 
supposed to contain only one particle of mass m. We 
aim at calculating formally the mean values (0(x, 0))bry 
(due to the translational y-invariance, y is set to for 
simplicity) of some, as yet unspecified, local field operator 
O. As is explained in section III.B, in the 1 + 1 Minkowski 
x-time and y-space, the boundary condition at a; = acts 
as the initial-time condition described by the boundary 
state 

\B) exp ||a+(0) + ^K{e)A+{-e)A+{9)^ |0) 

(5.1) 

which belongs to the bulk Hilbert space Ti. Thus, 



(O(x,0))bry = 



(O|0(x,O)|i?) 
(0|i?) ^ 



(5.2) 



note that the normalization is {Q\B) — 1. 

A systematic expansion for the one-point function l|5.2() 
can be obtained by using a complete system of the bulk 
n-particles states {\n)} forming TL as follows 

oo 

{Q\0{xMB) = Y^{Q\0{xMn){n\B). (5.3) 

ra=0 

In the rapidity representation, this formula reads 
(0|O(x,0)|B) = (5.4) 



(0|O(a;,0)|0> 



2ti 



2tt 



c\f) 

{o\o{x,omm) 

{o\o{x,o)\ei,02){ei,02\B) 



where \e) ^ A+ {0)\O) , ^2) = A+(^?i)A+(02)|O), etc. 

The matrix elements (0\B), (^i, 6*21-6), etc., which de- 
pend on the given boundary condition, can be calculated 
by using the explicit form of the boundary state H5.1|l . 
The normalization condition for the two-particle scatter- 
ing 

A{ei)A+{92) = 2^5{9i - 02) + S{02 - 0i)A+{02)A{0i) 

(5.5) 

implies that 

{0\B) ^ {Q\A{0)B) = |27r<5(0). (5.6) 



Analogously, we have 
{0i,02\B) = {Q\A{02)A{0i)B) 

poo -in 

= —K{0){2n)'S{0,+0)S{02-0), 



(5.7) 



and so on. 

The matrix elements (010(2;, 0)|6'), (010(2;, 0)|6li, 6I2), 
etc., which do not depend on the given boundary condi- 
tion, are known as the bulk multi-particle form factors. 



Their x-dependence can be factorized out by means of a 
translation on the operator O{x,0) plf: 

(O|O(x,O)|0) ^ e-^^^^'^^Fi, (5.8a) 

{O\O{x,O)\0i,02) = e-^<''°^'^<^^+''°^^<^^)F2{0i-02). 

(5.8b) 

and so on. The bulk form factors Fi = (0|O(0, 0)16*), 
F2{0i - 02) = (010(0,0)1611,612), etc. can be obtained 
explicitly in an axiomatic way, similarly to the case of 
the two-particle scattering ^-matrix. 

Finally, using the formulas H5.6|l - H5.8|l in the expan- 
sion (|5.4|) , the form-factor representation of the one-point 
function reads 



(O(2;,0))bry = (0|O(x,0)|0) + l^^ie-™^ 



/•oo 



I, ^^W^2(-2^)e^ 



(5.9) 



This formula is particularly useful for large distances x 
from the wall since it provides a systematic large-distance 
expansion. The first term is nothing but the bulk ex- 
pectation value of the operator O, (O), which is indeed 
dominant in the limit 2; ^ 00. The leading correction 
term comes from the one-particle state, with the particle 
mass m playing the role of the inverse correlation length. 
The next-to-leading correction term decays at large 2; as 
exp(— 2TOa;) multiplied by an inverse-power of mx^ and 
so on. 

The extension of the formula (|5.9() to general inte- 
grable 2D field theories with many-particle spectrum {a} 
is straightforward. We only write down the final result: 



(o(x,o))bry = (o) + ^|f(''V'- 

a 

+ 1 ^ K''\e)F!t'\-20)e 



(5.10) 



— 2raaX cosh t 



(a,b) 
rria—mb 



B. The electrical double layer 

Since the number (|3.18() and the charge (|3.19ll density 
profiles are determined by the boundary mean values of 
the exponential field, (e''?''"^("^'''))bry (q = ±1), the opera- 
tor of interest is 



Og(2;,0) = exp{iq&(?!)(2;,0)}, 



±1. 



(5.11) 



In the spectrum of the 2D bulk sine-Gordon theory, 
the first two lightest neutral particles are the Bi- and -82- 
breathers with the corresponding masses (see Eq. ()2.11|) ') 



mi 



m2 



2M sin 



\ 2 



2Msin(7rp), p< 



p<l, 
1 



(5.12a) 
(5.12b) 
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The mass m2 < 2toi, the equahty taking place in the 
DH hmit p — > 0. Consequently, when applying the large- 
X expansion l|5.1()|l to the boundary sine-Gordon theory 
with Dirichlct boundary conditions (|3.15|l . we have to 
consider the following leading terms 



(O,(x,0))bry = (O> 

- r^X"m("'(-2^;,)e 
Jo 



(5.13) 



+|F(^)(g)e-'"- + |Ff)(,)e- 



-2m\x cosh 6 



Here, the boundary couplings gi and §2 are given by 
Eqs. (|3.25al) and (|3.25bl) . respectively, and K'^^{9) = 



where the reflection amplitude R^^^ is given 



VT 

by Eq. 

Our preliminary task is to write down the bulk form 
factors in Eq. (|5.13(l for the exponential operator H5.11|l . 
We start with the Bi form factors which have been ob- 
tained in Ref. 22» The one-particle Bi form factor reads 



Fi\q) = -if?A(c'«'"^), q = ±1, 
where the parameter A is defined by 

1/2 / fP'' dt t 



(5.14) 



A 



2c„.(f) 



2 sin 



V 2 



cxp 



The two i?i-breathers form factor reads 



27r sin t 

(5.15) 



i"H(?;g) = -A2i?(0)(e'«^*), 



q = ±l, 



(5.16) 



where the function R{9) is given on the interval — 27r 
pTT < lm(0) < —piT by the integral 



i?(6l) = AAexp*! 8 
X sinli 
TV = exp <^ 4 



dt sinh(<) sinh(pt) sinh[(l + p)t\ 
t sinh^(2i) 

(5.17a) 

dt sinh(i) sinh(pt) sinh[(l -|- p)t\ 
1 sinh2(2t) 

(5.17b) 



t{l-- 



This function satisfies a useful relation 

sinh(6l) 



R{e)R{e ± ivr) 



sinh(6') ^ isin(p7r) ' 



(5.18) 



which, together with R{—0) = Sii{9)R{6), enables one to 
extend the definition of R{0) to arbitrary values of Im(6'). 
(Here 6*11 denotes the BiBi scattering matrix, see be- 
low). The evaluation of the one-particle B2 form factor 
can be based on a bootstrap procedure [21|. Namely, 
since the i32-breather is a boundstate of the two Bi- 
breathers (i.e. the BiBi scattering matrix has the B2- 
pole), the one-particle B2 form factor can be calculated 



from the two-particle Bi form factor (|5.16|) as follows |^ 



rFf^(g) = -ires,=oJ^; 



£=0^2 



(11) 



ipn - 



ipn 



2 ' 2 

(5.19) 

where F is related to the residue of the _B2-pole in the 
BiBi scattering as follows 



-i Tese=ipTrSii{9) = T^. 
Explicitly, one has 

smh{9) + isin(p7r) 



(5.20) 



^11 (0) = 



sinh(0) — isin(p7r) ' 
The one-particle B2 form factor thus reads 



r = v/2tan(p7r). (5.21) 



Ff^(q) = -A^ 



tan(p7r) 



1/2 



1 



R[m{l+p)] 



,iqb<f>\ 



(5.22) 



With regard to the bulk neutrality condition (e''"^) — 
(e^'*"^), under the transformation q —q the form factor 
F^'^\q) changes the sign, F^'^\+l) = -F^^\-l), while 



the form factors '{9;q) ((O^ and F^ ' (q) (f^:^ are 
unchanged, F^^^\9;+l) = F^^^\9;-l) and i^f ^(+1) = 



(2) 

Fl (—1). Consequently, as is clear from the definitions 
(|3.18l) and (|3.19|) . within the series representation H5.13|l . 
the term with F^*-^-* contributes to the charge density p 



while the terms with ^2^^'' and F^'^' contribute to the 
particle number density n. 



(2) 



C. Asymptotic charge profile 

The asymptotic large- a; behavior of the charge density 
is obtained in the form 



p{x) 



2n\ cos 



X exp 



(?) 



1 + cos 

dt t 
27r sin t 



(?) 

tanh(p(^) 



pTr\ 



X exp {—mix) , 



(5.23) 



where p = /3/(4 — /?) and the mass nii of the lightest 
i3i-breather is given by the relations H2.19|l and H2.20|l . 
The formula H5.23|l applies to the whole stability region 
f3 < 2 which is simultaneously the region of the exis- 
tence of the Bi- breather in the particle spectrum of 
the sine-Gordon model. With respect to the Poisson 
Eq. 1)3. 6|) . the deviation of the induced electrostatic po- 
tential from its bulk value H3.7|l behaves at large x as 



Sip{x) 



-^.-K p(x) I rn\. Since the pure exponential 



asymptotic decay of the DH results (|4.4|l is recovered, 
the concept of renormalized charge is applicable to the 
present model. The renormalized image charge tXren, de- 
fined by Eq. (|4.18|) with the replacement k — > mi, is 
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given by 



pn 



77/3 \l 2sin(]37r/2) 
dt t 

X exp 



1 + cos 1 



p7r\ 

2 y 



27r sin t 



tanh(p(/?). 



(5.24) 



It can be readily verified that the /3-expansion (|4.19|) is 
reproduced by this formula. In the limits /3 — > and 
if oo such that the product /3(p is finite, Eq. (|5.24|) 
reduces to the result H4.27(l of the nonlinear PB theory 
which is therefore adequate in such regime. For a given 
(3, increasing the potential difference ip to infinity, (Tron 
saturates monotonically at the finite value 



pTT 



77/3 \l 2sin(]37r/2) 
f"^ dt t 

X exp 



PTT\ . fpTT\ 
1 + COS I —:r I — sm I I 



2 / 



/Q 27rsint^ 
in agreement with the saturation hypothesis. 



V 2 /. 
(5.25) 



D. Asymptotic number density profile 

The asymptotic large- 2; behavior of the particle num- 
ber density n{x) is a more complicated topic. For the 
density deviation from its bulk value Sn{x) — n{x) — n, 
the form-factor asymptotic expansion (|5.13() gives 



Sn{x) - Sn'-^\x)+Sn^^^^x), 

X — >oo 

where the i?2-breather term reads 



(5.26) 



tan(p7r) 



1/2 



X exp(— m2x) 



i?[i7r(l+p)] 

(5.27) 



with 52 given by Eq. I)3.25b|l taken at 77 = i2i^, and the 
B\B\ integral term reads 

(5n(")(x) ^ / _;^(6i)e-2'"i^™''he^ (528a) 

x{e) = -]^n\^K^\e)R{e). (5.28b) 

The presence of the B2 term is restricted to p < 1/2 
(/? < 4/3), the BiBi is present in the whole stability 
region p < 1 {(3 <2). As soon as (3 (or, equivalently, p) 
has a strictly nonzero value, it can be shown from (|5.28(l 
that 



5n^^^\x) oc 



-2mix 



^mix 



at large x, 



(5.29) 



with 2mi > m2 for any /? > 0. The BiBi term thus 
becomes subleading in (|5.26(l . while the B2 term (15 .2711 
dominates: 



Sn{x) ~ Sn^'^^x), P > 0. 



(5.30) 



As concerns the /3 — > limit, considering i?(i7r) —1, 
A ^ 2^, .92 ^ Sv/V^I"^' + (^Vl6)], p ^ /3/4 and 
TO2 2k in H5.27|l leads to the expression 

Sn'-^^x) = P^n ^(^2 + exp(-2Kx) + 0{(3^), (5.31) 

which is twice larger than the expected result (|4.22() of 
the systematic /3-expansion. This means that the for- 
mula H5.3U|I does not reflect adequately the expansion of 
the asymptotic 5n{x) around the (3 = Q point. The rea- 
son for this inconsistency consists in the fact that, in the 
limit p ^ 0, besides the important equality of inverse 
correlation lengths m2 = 2mi = 2k, the value of x(^) ^ 
well as of all its derivatives go to infinity at 6* = 0. As a 
consequence, also the -Bi-Bi term (|5.28a(l . which is sub- 
leading for a strictly nonzero (3, contributes to the lead- 
ing order of the pure-exponential large- a; behavior [see 
the relation H5.35|l below]. It is important to add that if 
we would be able to perform all /3-orders of the large-x 
decay of S'^^^^x), we should arrive at the exact asymp- 
totic behavior (|5.29|) . valid for /? > 0, which is no longer 
purely exponential. This mathematical technicality was 
first observed in the study of finite-size effects for the 
(l-l-l)-dimensional sine-Gordon theory defined on a strip 
with Dirichlet-type boundary conditions .36.] : the next 
analysis follows a regularization procedure presented in 
Ref. [S^. The dangerous singularity of x(^), at 6* = in 
the limit p — » 0, can be isolated from x(0) in the following 
way 

cosh6' -I- cos(p7r/2) 

xiO) = —T-. tH^Xo ^ , 5.32 

coshe' — cos(p7r/2) 

where Xo(^) is a regular function of 9 around 6 = Q: 



Xo(^^)=Xo(0) + Xo(0)- 



(5.33) 



Here, Xo(0) corresponds to the classical treatment, Xo (0) 
to the first quantum correction, etc. Within the regular- 
ized form (|5.32l) complemented by the regular expansion 
(|5.33|) . the evaluation of the BiBi integral (|5.28a|) can be 
carried out in close analogy with Ref. l3(l In particular, 
for very large x and small p, one uses the results of [s^ 

d^ cosh6'-|-cos(p7r/2) ^_.^^^ ^^.^ e 
TV cosh 6* — cos(p7r/2) 

x-*ca sm(p7r/4) 

Similarly, one can derive that 

d6> cosh 6 + cos(p7r/2) 6>^ ^_2mx cosh e 
TT cosh0 — cos(p7r/2) 2 



~ 2 sin ( - — 1 cos 

X — ''OO V 4 / 



r^Vn«(^)e-^-. (5.34b) 



Regarding the explicit forms of Xo(0) and Xo(0)j the re- 
sult for the BiBi integral is such that 



(5n(")(a:) 



lim 



(5.35) 
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In view of relation (|5.31(l . the sum in Eq. (|5.26() thus 
reproduces the needed result (|4.22l) . To conclude, the 
small /3-expansion around the /3 = point of the asymp- 
totic number density profile involves artificial contribu- 
tions from the BiBi integral term, this term being to- 
tally absent at a strictly nonzero (possibly very small) 
value of (3. In other words, the small-/? expansion of the 
asymptotic density profile, when truncated at some finite 
/3-order, does not reflect adequately the asymptotic den- 
sity profile at a nonzero value of /3. From this point of 
view, the DH results for number density profiles have to 
be taken with caution also for other boundary Coulomb 
systems. 



VI. THE FREE-FERMION /3 = 2 POINT 

For the sake of completeness, we summarize in view of 
the subjects of present interest the exact results for the 
2D electrical double layer at the free-fermion f3 = 2 point 
of the Thirring representation of the Coulomb gas [T3| . 
The exact solution for the species density profiles at an 
arbitrary distance x > from the ideal-conductor wall 
reads 



m 

n±{x) - n± = — 

ZTT 



dl 



m 



Ki exp(±/3(^) -I- m 



m cosh(/3(^) 



X exp(— 2K/a;), 



(6.1) 



where m — 2ttz is the rescaled fugacity (equal to the 
soliton mass M, see formula (|2.17|l taken at = 1/2 
and p = 1), Ki — {m? + Py/"^ and n± = n/2 are the bulk 
densities regularized by considering a hard-core repulsion 
around each particle. The short-distance limit of Eq. 
(EH), 



n+[x] ^ — as X 
^ ^ 2x 



0, 



(6.2) 



is of the expected form H3.12|l . The charge density, cal- 
culated from Eq. (|6.1|) . implies via the Poisson equation 
the following deviation of the electrostatic potential from 
its bulk value: 



X f \ Ai ^ msinh(2(^) 

Sifix) = - / dl- exp{-2Kix) 

Jq 2ki mcosh^Zip) + Ki 

At asymptotically large distances from the wall. 



(6.3) 



S(p{x) 



1 / TT \l/2 

-tanhiysl exp(— 2r7ia;). (6.4) 

4 \mxJ 



This asymptotic behavior differs fundamentally from the 
purely exponential DH prediction 14.4a|l . The concept of 
renormalized charge is therefore not applicable at the free 
fermion point. The reason for the fundamental difference 
is obvious. The lightest i?i-breather disappears from the 
particle spectrum of the sine-Gordon model just at the 
free-fermion point (3 = 2, and the asymptotic behavior 
of 5lp{x) starts to be governed by the soliton-antisoliton 



pair. Since mi 2M (= 2m) as /3 — > 2, the par- 
ticle mass in the exponential decay on the rhs of Eq. 
(|6.4|) is a continuous function of /? at = 2. On the 
other hand, the position-dependent prefactor [mx)~^/'^ 
in the formula (|6.4() , determined by the form- factor of the 
soliton-antisoliton pair, has no "continuous" analogue in 
the leading asymptotic behavior of 5lp{x) inside the sta- 
bility region /3 < 2. The basic qualitative features of the 
results at the [3 = 2 point are expected to be present 
also for /? > 2, up to the Kosterlitz-Thouless point /3 = 4 
where the 2D sine-Gordon theory ceases to be massive. 
In the limit ip oo, Sip{x) of Eq. (|6.3I) saturates at 



S(p*{x) 



1 



'-Ko{2m,x) 



(6.5) 



This function is finite in the whole electrolyte region x > 
0, which confirms the validity of the hypothesis of the 
electric potential saturation |10| . 



VII. CONCLUSION 

The main aim of this paper was to test basic concepts 
used in the theory of highly asymmetric Coulomb flu- 
ids on the exact solution of a 2D electrical double layer. 
This model is mappablc onto the 2D semi-infinite sine- 
Gordon field theory with Dirichlet boundary conditions 
which do not break the intcgrability property of the bulk 
sine-Gordon model with the known particle spectrum. At 
large distances from model's interface, the induced elec- 
tric potential has the pure exponential decay for small 
enough inverse temperatures (couplings) (3 < 2, includ- 
ing the DH — !■ limit. This fact confirms the ad- 
equacy of the concept of renormalized charge for weak 
couplings. In the extreme case of an infinite potential 
difference between model's interface and the bulk inte- 
rior of the electrolyte, the renormalized charge saturates 
at a finite value which is in agreement with the saturation 
hypothesis. 

As concerns the applicability of the present results in 
the context of "real" 3D Coulomb fluids, we expect the 
reduced dimensionality to be unimportant because, as 
soon as one uses the 2D harmonic (i.e., logarithmic) po- 
tential, basic screening sum rules for the charge-density 
averages are maintained. Another "imperfection" of our 
model is the absence of a hard core around each charged 
particle: in 3D, the hard-core repulsion is inevitable to 
avoid the thermodynamic collapse at any (3 > 0. Since 
the charge sum rules arc not affected by a short-ranged 
regularization of the Coulomb potential 12j, also this 
model's simplification seems to be irrelevant. 
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